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^ HOW TO USE MOLECULAR CLOUDS TO STUDY THE 

O PROPAGATION OF COSMIC RAYS IN THE GALAXY 

Observations of molecular clouds in the gamma ray domain provide us with 
a tool to study the distribution of cosmic rays in the Galaxy. This is because 
cosmic rays can penetrate molecular clouds, undergo hadronic interactions in 
the dense gas, and produce neutral pions that in turn decay into gamma rays. 
^ The detection of this radiation allows us to estimate the spectrum and intensity 

^ of cosmic rays at the cloud's position. Remarkably, this fact can be used to 

constrain the cosmic ray diffusion coefficient at specific locations in the Galaxy. 
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^ Galactic Cosmic Rays (CRs) arc believed to be accelerated by SuperNova 

ON Remnants (SNRs).i The main (but not conclusive) arguments supporting 

this idea are the fact that SNRs can provide the total power required to 
maintain the galactic CR poptilation, and the belief that an efficient accel- 
' eration mechanism, diffusive shock acceleration,^ operates in these objects. 

The detection of some SNRs in TeV gamma rays,^ expected if they indeed 
are the sources of CRs,* also supports this scenario, though the emission 
\ l might have a leptonic origin and not be related to the acceleration of CRs.''' 

^ Besides the total CR power of the Galaxy, a theory of CR origin must 

also explain the spectrum, isotropy, and chemical composition of CRs.^ This 
requires a thorough knowledge not only of the nature of CR sources, but 
^ also of the way in which CRs propagate in the Galaxy. This is because the 

observed properties of CRs results from the equilibrium between the injec- 
tion rate of CRs from sources and their diffusive escape from the Galaxy. 

Measurements of the amount of spallation suffered by CRs allow us to 
infer the average residence time of a CR of energy Ep in the Galaxy as 
tres oc E~^, with 6 ~ 0.3—0.7. If h is the length a CR has to move away from 
its source before escaping the Galaxy (i.e. the Galaxy's thickness), then the 
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diffusion coefficient reads: Dgai ~ h^^ /Ues ~ lQ'^^{Ep/l{) GcYf cm^/s.. 
However, this has to be intended as the average diffusion coefficient in the 
Galaxy, and local variations (both in time and space) might exist. 

In particular, the diffusion coefficient might be suppressed close to CR 
sources. This is because CRs can excite magnetic turbulence while stream- 
ing away from their acceleration site. This would enhance the scattering 
rate of CR themselves and thus reduce the diffusion coefficient.^ The prob- 
lem of estimating, on theoretical grounds, the diffusion coefficient around 
CR sources is far from being solved, mainly because of its intrinsic non- 
linearity and because various mechanisms might damp the CR-generated 
waves and thus affect the way in which CRs diffuse. ^'^ 

In principle, gamma ray observations can provide us with constraints on 
the diffusion coefficient close to CR sources. Once escaped from their 
sources, CRs undergo hadronic interactions with the surrounding gas and 
produce gamma rays. The cliara(;teristics of such radiation (in particular its 
spectrum and intensity as a function of the time elapsed since CRs escaped 
the source) depend on the value of the diffusion coefficient that can thus be 
constrained, if a reliable model for CR acceleration at the source is avail- 
able. The presence of massive Molecular Clouds (MCs) close to the source 
would enhance the gamma ray emission, making its detection more proba- 
ble. Studying such radiation is of great importance not only in order to reach 
a better understanding of how CRs diffuse, but also bc;eause its detection 
can provide an indirect way to identify the sources of galactic CRs.^°"^^ 

Molecular clouds as cosmic ray barometers 

Consider a MC with mass M^i at a distance d, located in a region in the 
Galaxy where the CR (protons) intensity is Jcr- For illustrative purpose, 
we set here Jcr — KcrE~°'. We further assume that the high energy CRs 
(the ones with energy above the threshold for 7r°-production) can freely 
penetrate the cloud. Under these assumptions the expected gamma ray 
flux from the MC due to proton-proton interactions is given by: 



where app ~ 34 mb is the interaction cross section, nip is the proton mass, 
and Yy depends on a and is tabulated in Ref. 6. Equation 1 is valid at high 
energies only {E^ > 10 GeV), while at lower energies the spectrum (in log- 
log scale) is symmetric with respect to the energy E^ = m^a /2 ^ 70 MeV. 
Assume now that the CR intensity in the region under exam differs by 
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a factor 5{Ep) = JcR{Ep) / Ji,g{Ep) from the one measured at the Earth 
(which is:!^ {'i'K/c)E^Jbg{Ep) ~ 6 x 10-3(i;p/TeV)-°-^eV cm'^), so that 
Equation 1 can be rewritten as: 

E^^F,{E,) ~ 2.5 X 10-" C^jI^) (^^) S{Ep)E-^-'TeV/cmys (2) 

where M5 is the mass of the MC in units of lO^M©, dkpc is the distance in 
kpc, and /o.i = (//O.l) « 1 takes into account the fact that on average CRs 
with energy Ep produce gamma rays with energy E^ f x Ep. Moreover, a 
multipHcative factor of 1.5 has been appHed to account for the contribution 
to the emission from nuclei heavier than H both in CRs and in tlic MC gas. 

If gamma rays are detected from a MC, and its mass and distance are 
known from its CO emission, then Equation 2 allows one to measure 
both the spectrum and intensity of CRs at the MC's location, and thus 
the quantity S{Ep). It follows that MCs can be effectively used as probes 
of the energy density of CRs at different locations in the Galaxy, and for 
this reason have been sometimes referred to as CR barometers Since 
observations in the GeV range of the Galaxy suggest that, on large spatial 
scales, CR variations are not very large, ^^"^'^ ^(Ep) is normally interpreted 
as the excess above the galactic CR background. Some examples of how 
to use this fact to estimate the CR diffusion coefficient are given in the 
following sections. 

The galactic centre ridge 

The H.E.S.S. collaboration reported on the detection of diffuse gamma ray 
emission from the Galactic Centre (GC) ridge. Remarkably, the emission 
correlates spatially with a complex of giant MCs, suggesting that the emis- 
sion itself is likely to be originated by CR hadronic interactions in the dense 
ambient gas. The TeV gamma ray spectrum extracted from the region of 
galactic coordinates \l\ < 0.8°, |6| < 0.3° can be fitted by a power law 
with index a ^ 2.29 and normalization ^ 1.73 x IQ-^TeV^^cm^^s^^sr^^. 
At the distance of the GC (~ 8.5 kpc) the emitting region has a size of 
w 240 X 90 pc, and encloses a total gas mass of 1.7 - 4.4 x lO'^ Mq.^^ 

By assuming that all the observed gamma ray emission has an hadronic 
origin, we can use Equation 2 to show that, close to the GC, the CR intensity 
at 10 TeV is enhanced with respect to the local one by a factor of (5 ~ 3 — 9, 
and grows at higher energies due to the hardness of the measured emission. 
Moreover, the total energy in form of CRs needed to explain the emission is 
« lO^^erg, if the observed spectrum is extrapolated from 1 GeV to 1 PeV. 
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Thus, a single SNR with explosion energy Esn ~ lO^-'^erg and acceleration 
efficiency 77 « 10% might explain the whole CR excess in the region. 

One possible CR accelerator is the SNR Sgr A East, located close to the 
GC and with an estimated (but uncertain) age of tage ~ 10 kyr. If the CRs 
responsible for the observed emission are assumed to be released during the 
earliest phase of the SNR evolution, then it is possible to estimate the 
diffusion coefficient by requiring them to ffil the emitting region of size 
w 120 pc in a time tage- This gives D « tage ~ lO^^cm^/s, which, 
taking into account the large uncertainties in the assumptions made, is in 
rough agreement with (maybe somewhat smaller than^*) the value of the 
diffusion coefficient of ^ TeV CRs in the galactic disk. The presence of a 
massive MC at I = 1.3° (^ 200 pc) which does not show prominent gamma 
ray emission suggests that CRs accelerated at the GC have not diffused yet 
up to that distance, and this justifies our assumption for Z-y.^^ 

This estimate of the diffusion coefficient has to be regarded as a very 
rough one, due to the presence of many uncertainties in its derivation. For 
example, another source (e.g. the black hole Sgr A*^^), or many sources (e.g. 
an ensemble of SNRs^^) might be the accelerators of the CRs responsible 
for the TeV emission. Moreover, proton energy losses (ignored here) might 
affect the estimate of the diffusion coefficient by a factor of « 2.-^^ 

Better constraints on the diffusion coefficient and on the nature of the 
source(s) of the accelerated particles could be obtained if the spatial vari- 
ation of the spectrum of the gamma ray emission was known. Measuring 
such variation, though challenging for currently operating instruments, is 
within the capabilities of future ones like the Cherenkov Telescope Array. It 
seems thus appropriate to discuss the future perspectives in this direction. 

According to the scenario summarized above, the SNR Sgr A East, 
located close to the GC, released about 10^*^ erg in form of CRs about 
10 kyr ago. Since high energy CRs are believed to be accelerated by young 
SNRs (aged up to few thousands years), we are allowed to consider this as 
an impulsive injection event. If CRs are injected with a power law spectrum 
DC E~'^, their time dependent spatial distribution around the source is:^^ 

WE,,t) = ^3/.i,(^4r/^^..„) ^ e-i^S (3) 

where Ep^rnax{Ep,min) is the maximum (minimum) CR energy, = 
■^/4:D{Ep) t is a diffusion length and D the diffusion coefficient. The CR 

intensity is the sum of Equation 3 and the galactic CR background. 

The expected excess over the background for CRs with energy 3 TeV is 
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Fig. 1. Upper panel: excess of 3 TeV CRs over the CR background. Lower panel: slope 
of the gamma ray cmissivity at 300 GcV for S = 0.3 (solid) and 0.4 (dashed). Horizontal 
lines represent the observed slope (measured value it error). Dotted vertical lines delimit 
the emitting region. Arrows refer to the positions of the peaks in the gas distribution. 

shown in Figure 1 (top panel) as a function of galactic longitude. These 
arc the CRs that produce most of the emission detected above ^ 300 GeV 
in the region \l\ < 0.8° (marked by dotted vertical lines in Figure 1). The 
slope of the gamma ray emissivity at 300 GeV and for a diffusion coefficient 
D oc with 6 = 0.3 is plotted as a solid line in the bottom panel, where it is 
compared with the measured one, represented by the horizontal dotted lines 
(measured slope ± statistical error). The arrows represent the (projected) 
position of the four prominent peaks in the gas density distribution. If the 
total gas mass is concentrated at these locations and if projection effects 
in the gas distribution do not play an important role, the spectral slope of 
the emissivity should roughly represent the observed slope also. 

The results in Figure 1 can be interpreted as follows. Within roughly one 
diffusion length, the gamma ray spectrum hardens with the distance from 
the source. This is because higher energy CRs occupy a larger volume of 
space around the source, and thus the low energy emission is progressively 
suppressed at larger and larger distances. On the other hand, at distances 
much larger than one diffusion length, the CR excess vanishes and the steep 
spectrum characteristic of the CR background is recovered. Though these 
calculations are far too crude, they might still serve to estimate the order of 
magnitude of the expected spectral variations within the emitting region. 
From Figure 1 one sees that variations of the order of Aa ~ 0.1 — 0.2 are ex- 
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pcctcd for a diffusion coefficient with slope 5 = 0.3. A stronger dependency 
on energy would produce more pronounced variations, as illustrated by the 
dashed line which refers to ^ = 0.4. Values of 6 significantly larger than 0.3 
would require a harder injection spectrum of CRs in order to preserve the 
observed slope in gamma rays. Such variations in the slope of the emission, 
if detected, would support the scenario of a localized injection of CRs close 
to the GC and provide better constraints on the diffusion coefficient. Com- 
plementary studies of the three dimensional distribution of the gas in the 
region are also needed in order to control projection effects. 

Molecular clouds in the W28 region 

W28 is a SNR in its radiative phase of evolution, located at a distance of 
~ 2 kpc, in a region rich of dense molecular gas. Gamma ray emission has 
been detected from the surroundings of W28 both at TeV^® and GeV en- 
gj.gies,30.3i by HESS, FERMI, and AGILE, respectively. The TeV emission 
correlates quite well with the position of three massive (« lO^M©) MCs, 
one of which is interacting with the north-eastern part of the shell (and 
corresponds to the TeV source HESS J1801-233), and the other two being 
located to the south of the SNR (TeV sources HESS J1800-240 A and B) . 

Let us assume that this gamma ray emission is the result of hadronic 
interactions of CRs that have been accelerated at the SNR and then escaped 
in the surrounding medium. "^^ In this scenario, the distribution of runaway 
CRs around the SNR is still given by Equation 3. A diffusion coefficient 
D = X-Dgai E^'^ is assumed, where x represents possible deviations 
with respect to the average diffusion coefficient in the Galaxy. Equation 3 
tells us that, up to a distance equal to the diffusion radius Ra, the spatial 
distribution of CRs around the source is roughly constant, and given by 
^CR rjEsN I R% On the other hand, the observed gamma ray flux from 
each one of the MCs is: oc NcRMd/d^ ■ By using the definitions of 
NcR and Rd we can finally write the approximate equation, valid within a 
distance Rd from the SNR: 



Estimates can be obtained from observations for all the physical quantities 
in Equation 4 except for the CR acceleration efficiency 77 and the diffusion 
coefficient x^gai- By fitting the TeV data we can thus attempt to constrain 

a combination of these two parameters (namely ij/x^^^)- Given all the un- 
certainties above, our results have to be interpreted as a proof of concept 
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Fig. 2. Simultaneous fit to the three TeV sources detected by HESS in the W28 region. 
Figure from Ref. 32. See text for more details. 

of the fact that gamma ray observations can serve as tools to estimate the 
diffusion coefScient. More detections of SNR/MC associations are needed 
in order to check whether the scenario described here apphcs to a whole 
class of objects and not only to a test-case as W28. Future observations 
from the Cherenkov Telescope Array will most likely solve this issue. 

Figure 2 shows a simultaneous fit to the HESS data for the three massive 
MCs in the W28 region, obtained by setting ~ 20. This implies that 

the normalized diffusion coefScient x has to be much smaller than 1 for any 
reasonable value of 77 < 1. For example, an acceleration efficiency 77 = 30% 
corresponds to x = 0.06, which in turn gives a diffusion distance for TeV 
particles of Rd « 60 pc. This means that the results in Figure 2 are valid 
if the physical (not projected) distances between the MCs and the SNRs 
do not significantly exceed Rd- Small values of the diffusion coefficient have 
been also proposed in Ref. 31,33,34. Note that, since we are considering 
gamma rays in a quite narrow energy band around ~ 1 TeV, we are actually 
constraining the diffusion coefficient of CRs with energy Ri 10 TeV. With 
some caveats, "^^'^^ observations in the GeV range might be used to constrain 
the diffusion coefficient down to GeV particle energies. 

IC 443 and other examples 

The SNR IC443 has an estimated age of tage ~ 3 x 10^ yr, explosion energy 

of EsN ^ 4 X 10'''°crg, and distance « 1.5 kpc.'^^ It has been detected at 
both GeV^'' and TeV"^* energies. A power law fit to the TeV emission is: 
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- 10-11 (£;^/0.4TeV)-3-icni-2s-iTcV-^ Observations of the CO line 
indicate that about 10^ of molecular gas are present in the region. 
Though the GeV emission might be explained as the emission from shocked 
clouds,^^ this scenario cannot account for the very high energy CRs that 
originate the TcV photons. If these photons are produced by runaway CRs 
we can attempt to estimate the diffusion coefficient, as done for W28. If the 
TeV radiation is produced in a MC of mass ~ 10^ , then from Equation 2 
wc can infer an excess in the intensty of multi-TeV CRs of the order of 
~ 100. If these particles had time to diffuse up to a distance Rd, their 
total (integrated over volume) energy is ■K^^'^EpNcR{0, Ep,tage)Rd ~ 3 x 
10*^(iid/10 pc)^ erg (see Equation 3). This number has to be multiplied 
by « ln(ii'™"^/£^™*") to account for the extension of the spectrum of CRs 
injected by the SNR and then equated to rj Esn- This gives (very roughly) 
a diffusion coefficient for multi-TeV CRs of D w lO^^rj/O.l)^^ cm^/s. 
This is, again, much smaller than the average diffusion coefficient in the 
Galaxy,^^ and corresponds to Ra ~ 20 pc. 

Other SNR/MC associations detected in TeV gamma rays that would 
deserve attention include CTB 37A, HESS J1745-303, and W51.^° The 
approach described here might now be extended to the GcV range, where 
several SNR/MC associations are being detected by FERMI. The Cygnus 
region, rich of dense gas and detected at multi-TeV energies, also promises 
to be an ideal region of the sky for this kind of studies. 



Conclusions 

MCs, when detected in gamma rays, can serve as probes of the CR intensity 
throughout the Galaxy. If they are located close to CR sources, the intensity 
of the gamma ray emission might impose constraints on the local diffusion 
coefficient. In at least one case (SNR W28), preliminary evidence has been 
reported for a significant suppression of the diffusion coefficient with respect 
to the average galactic one. This suppression might be the result of an 
enhancement in the magnetic turbulence due to the streaming of CRs away 
from the source. Observations with future facilities such as the Cherenkov 
Telescope Array will conclusively test the feasibility and reliability of the 
approach to constrain CR diffusion that has been presented here. 
Support from the EU is acknowledged [FP7 - grant agr. n° 256464]- 
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